Noisy NF-ftB oscillations stabilize and sensitize cytokine signaling in space. 
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NF-ftB is a major transcription factor mediating inflammatory response. In response to pro- 
inflammatory stimulus, it exhibits characteristic response - a pulse followed by noisy oscillations 
in concentrations of considerably smaller amplitude. NF-^B is an important mediator of cellular 
communication, as it is both activated by and upregulates production of cytokines, signals used by 
white blood cells to find the source of inflammation. While the oscillatory dynamics of NF-^B has 
been extensively investigated both experimentally and theoretically, the role of the noise and the 
lower secondary amplitude has not been addressed. 

We use a cellular automaton model to address these issues in the context of spatially distributed 
communicating cells. We find that noisy secondary oscillations stabilize concentric wave patterns, 
thus improving signal quality. Furthermore, both lower secondary amplitude as well as noise in the 
oscillation period might be working against chronic inflammation, the state of self-sustained and 
stimulus-independent excitations. 

Our findings suggest that the characteristic irregular secondary oscillations of lower amplitude are 
not accidental. On the contrary, they might have evolved to increase robustness of the inflammatory 
response and the system's ability to return to pre-stimulated state. 



INTRODUCTION 

The regulatory network of NF-^B is an important 
constituent of the immune system as it regulates hun- 
dreds of genes in response to extracellular stimuli such as 
pathogens, cytokines and stress [U[2]. These responses 
include apoptosis, cell proliferation, and inflammatory 
response [TJ [3] . 

During inflammatory response, the NF-^B network is 
activated by an increase in extracellular concentration 
of cytokines; small signaling molecules commonly used 
in intercellular communication. Once activated, NF-^B 
up-regulates the cells' own cytokine production, thereby 
amplifying the external signal and passing it to the neigh- 
bouring cells [4j [5] either by diffusion or through gap- 
junctions - channels formed by physically interacting 
cells [5]. 

When passing the signal from one cell to the next in 
this manner, the tissue acts as an excitable medium. Re- 
cent theoretical research has shown that propagating ele- 
vated cytokine concentrations as waves through this "ex- 
citable tissue" might represent an optimal way of passing 
the signal to the blood vessels, where the cytokines are 
absorbed [6j [7] . When neutrophils (white blood cells) de- 
tect the cytokines in the blood stream, they will follow 
the cytokine concentration gradient back to the site of 
infection [8], in order to destroy or contain the cause of 
the infection [9]. 

Single cell measurements [4] [10] and modeling ap- 
proaches [TTJ[T2] revealed that the concentration of active 
NF-ftB oscillates with an initial high amplitude peak and 
several consecutive lower amplitude peaks. The experi- 
ments also revealed a considerable amount of noise in 
this response, and simulations suggest that the noise is 
in fact induced by an inherent component of the network 



[T3HT6] . The effect of noise on the wave propagation of 
cytokines through the tissue has not yet been elucidated, 
but the presence of a noise-inducing component suggests 
that the noise itself may play a key role in the immune 
response. This has motivated us to ask the following 
questions: How does noise affect wave propagation? Is 
the system equally sensitive to irregularities in the pe- 
riod of oscillations as to irregularities in the refractory 
period? Does noise contribute to the onset of chronic 
inflammation - the state where waves are self-sustained 
and do not depend on stimulus? 

THE MODEL 

The tissue is modeled as an excitable media using 
a parallel cellular automaton, an algorithm frequently 
used as a mathematical idealization of biological self- 
organizing systems [17-21 j. Our model comprises 101 x 
101 cells placed on a square grid. All cells are initially 
inactive, but once activated by their neighbouring cells 
they start producing cytokines in an irregularly oscillat- 
ing fashion with an initial high peak followed by several 
lower peaks (see Fig. [IJi). The amplitude of the initial 
peak is normalized to 1 and the mean oscillation period 
is estimated from current literature to be 100 minutes 
[13], corresponding to 12 time steps in the model. The 
amplitude A of all secondary peaks and the relative size 
T] of the gaussian noise in the oscillation period are kept 
as free parameters, but retained similar to values exper- 
imentally measured in nuclear NF-^B. 

A cell is activated if the combined cytokine production 
of its 8 nearest neighbours, the Moore- neighbour hood, 
exceeds a threshold 6 = 3. That is, a cell can be ac- 
tivated if 3 of its neighbours are in their first peak of 
cytokine production, or if e.g. 2 neighbours are in their 
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FIG. 1. A: The oscillatory cytokine production of a single 
cell subject to constant stimulation. The amplitude of the 
initial peak is normalized to 1 while all secondary peaks have 
amplitude A. When the cell leaves refractory period, it is 
here immediately re-excited to its initial peak. B-E: show 
wave propagation in an ideal system without any noise. F-I: 
have r\ — 40 % noise in the oscillation period and no noise in 
the refractory period, r\ r — 0. J-M: show the system with no 
noise in the oscillation period, r\ Q — 0, but r\ r — 10 % noise in 
the refractory period. The secondary amplitude is A — 0.25 
in all panels A-M. In the time series F-I and J-M, snapshots 
are taken at t = 20, 60, 100 and 200 hours. We see that both 
kinds of noise gradually generate spirals of active cells, which 
may mislead neutrophils and cause chronic inflammation. 



first peak and 1/A other neighbours are in a secondary 
peak. When a cell is activated, it enters a refractory 
period in which it is completely unresponsive to the lo- 
cal cytokine concentration, and will produce cytokines at 
its own pace. The mean refractory period is estimated 
from current literature to be 200 minutes [22 j, but here 
is also added a gaussian noise of relative size r? r , kept 
as a free parameter. Once out of the refractory period, 
the cell becomes susceptible to input signals again, and 
if the local cytokine concentration should rise above the 
threshold, the cell will re-initiate the cytokine produc- 
tion, responding with the initial high peak followed by 
several low peaks (see Fig. [IJi). 

A stimulus of 3 x 3 constantly excited cells in the middle 
of the grid was added to act as a pathogen and initiate a 
train of cytokine waves propagating through the tissue. 



An example of the effect of noise in the wave propa- 
gation is shown in Fig. [I] Here, noise is added to the 
period of oscillations (Fig. JTt^-I) and to the duration of 
the refractory period (Fig. [llj-M). The addition of noise 
has two biologically interesting effects on the system: 



Higher secondary amplitude renders the system 
more robust to noise. 



In order to destroy or contain pathogens, the neu- 
trophils need a distinct concentration gradient to follow 
to the site of infection. Adding sufficient amount of noise 
to the period of oscillations and/or in the length of the re- 
fractory period of each cell causes the propagating waves 
to disintegrate or to lose their radial gradient. Biolog- 
ically, this may compromise the effectiveness of the in- 
flammatory response. 

The time it takes the system to lose its wave struc- 
ture depends on the amount of noise in the system. In 
order to examine this time before breakdown, the state 
of the system over time is compared to the state of the 
system the first time the waves fill the entire grid. This 
is achieved by calculating the Pearson correlation coeffi- 
cient between cytokine concentration of the first grid in 
which the waves have reached the outer edges and the 
grid in all later time steps: 



Ci,t = ( 



9i -(gi) 9t-(gt) 



(i) 



Here, gi is the cytokine levels of the first grid that is 
filled with waves, and g t is the grid at a later time step. 
Angle brackets denote spacial average and a denotes spa- 
tial standard deviation. The correlation coefficients are 
plotted as a function of time in Figj2]A.-B. The correla- 
tion coefficients decay due to noise destroying the waves. 
The characteristic time for breakdown of wave structure 
is defined as the last time the correlation coefficient to 
the first wave is more than 0.25 (see Fig. [2^-B). The 
exact value of this threshold is not essential for the main 
results, as long as it is low enough to allow for some ir- 
regularity in wave patterns. 

The correlation method allows us to investigate the 
characteristic time for breakdown as a function of noise 
in the period the oscillations and noise in the refractory 
period. Interestingly, the secondary amplitude A, has a 
significant effect on the systems ability to sustain wave 
propagation. As is seen in Fig. [2p-E, the parameter 
range with stable wave patterns increases with increas- 
ing secondary amplitude. Thus, high secondary ampli- 
tude makes the system more robust to noise. One pos- 
sible explanation for this is that secondary waves might 
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FIG. 2. A-B: Correlation to the first wave as a function 
of time in hours, illustrating how the characteristic time for 
breakdown of wave propagation is found. The red line cor- 
responds to the correlation threshold of 0.25. A shows cor- 
relations of a system with little noise that is able to sustain 
waves. Parameters: A = 0.25, r\ Q — 0.02% and r\ r — 0.01%. 
B shows correlations of a system with sufficient noise to break 
the waves. Parameters: A — 0.25, r\ Q — 0% and r\ r — 0.05%. 
C-E: Characteristic time for breakdown as function of noise 
levels in the period of oscillations and in the length of refrac- 
tory period for different secondary amplitudes. The secondary 
amplitude in panel C-E is 0.15, 0.25, and 0.35, respectively. 



help straighten up the lagging sections of the major wave- 
fronts by increasing the cytokine concentration at the in- 
flections points and effectively accelerating them. Since 
lagging cells can get excited by having e.g. two neigh- 
bours in primary peak and 1/A neighbours in secondary 
peak, a high secondary amplitude will lead to a larger 
acceleration. There is, however, an upper limit to the 
secondary amplitude: Waves should only originate from 
the central stimulus, so the combined secondary ampli- 
tudes from the 8 nearest neighbours of a cell should never 
exceed the activation threshold = 3. That is, if the am- 
plitude of secondary oscillations is larger than 0/8, cells 
have a possibility of getting activated by the secondary 
oscillations of its neighbours. This could happen any- 
where in the tissue, not just close to a stimulus, causing 
an inflammatory response from healthy tissue leading to 
chronic inflammation. Although this upper limit on the 
secondary amplitude is unknown for living cells, it is in- 
teresting to note that the experimentally observed sec- 
ondary amplitude is significantly smaller than the initial 
peak, with the secondary amplitude being about 20 % of 
the first peak [T3] ! 



Noise leads to self-sustained excitations resembling 
chronic inflammation. 

To further address the issue of chronic inflammation, 
we tried removing the stimulus after 120 periods of os- 
cillations, corresponding to 200 hours. If this eliminated 
the cytokine production, such that no cell in the tissue 
was re-excited to its first peak after another 48 periods of 
oscillations, we defined the system to be sensitive to the 
removal of the stimulus. If not, the system was declared 
insensitive. A system that is insensitive to the removal of 
the stimulus is not able to relax back to the pre-infected 
state, and is reminiscent of chronic inflammation. 

We found that noise in the oscillation period and/or in 
the duration of refractory period can create stable addi- 
tional sources and spirals, which excite the system inde- 
pendently of the stimulus, see Fig. [5JA.-D. Thus, noise not 
only destroys the radial gradient of the wave pattern, so 
that neutrophils will be unable to locate the stimulus, but 
the system can become insensitive so that removing the 
stimulus will not terminate the inflammatory response. 

As the system is particularly sensitive to noise in the 
refractory period, we investigate how the insensitivity de- 
pends on the level of refractory noise, r] r . We define the 
probability of becoming insensitive to the stimulus, 
as the fraction out of 30 runs, where the system becomes 
insensitive. 

Interestingly, with no noise in the oscillation pe- 
riod, the probability to become insensitive is a non- 
monotonous function of refractory noise, as seen in see 
Fig. [3j5. With no noise, the system is sensitive, P(I) = 
at T] r = 0, but a small increase in noise has a dramatic 
effect and completely destroys the sensitivity, P(I) ~ 1 
when n r = 0.01. At this small noise level, the insensitiv- 
ity is caused by a double spiral very close to the center, 
where the stimulus once was, as seen in Fig. [3JA.-B. To 
sustain this double spiral, the previous wavefront repeat- 
edly excites the central cells through a narrow passage 
of cells, as illustrated in Fig. [3^3. If the noise level is in- 
creased further to 0.04 < r\ r < 0.1, this signal can not be 
transmitted through the narrow passage, and the prob- 
ability that the system is insensitive decreases. This ex- 
plains the decline in probability for insensitivity in figure 
[3j3. Increasing noise further, n r > 0.1 produces spirals 
in rich numbers everywhere in the tissue, as shown in 
Fig. [3p-D. These spiral structures are quite stable and 
can continue to excite cells on long timescales after the 
stimulus has been removed. 

To better understand the non-monotonous behavior 
of the systems sensitivity, we have tested the case with 
larger system size, which should give higher chance for 
spirals to originate away from the center. As expected, 
Fig. [3^ shows that the onset of permanent insensitiv- 
ity has been expedited to lower noise levels, reducing the 
range of stimulus sensitivity. Furthermore, adding noise 
to the oscillation period prevents the formation of double 
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FIG. 3. A-D show two distinct ways the system can be- 
come insensitive to the stimulus. A-B: For small noise levels, 
a double spiral appears close to the removed stimulus. To 
sustain this double spiral, the previous wavefront needs to 
excite the central cells through a narrow passage. Parame- 
ters: A = 0.25, r\ = 0%, r\ r = 2%. C-D: For larger noise 
levels spirals arise numerous places in the tissue. Parameters: 
A = 0.25, J] = 0%, J] r = 20%. E: The probability P(I) that 
the noise in the refractory period causes the system to become 
insensitive to removal of stimulus as a function of noise in the 
refractory period. Interestingly, the curve is non-monotonic. 
F: Same as E but with a grid size of 202 x 202 cells instead 
of 101 x 101. The larger system size expedites the emergence 
of stable spirals. G: Same as E but with with a noise in the 
oscillation period of rj Q = 10 %, which destabilizes the dou- 
ble spirals for low noise levels, making the insensitivity curve 
monotonic. 



spirals close to the center, as the cytokine signal is again 
prevented to pass through the narrow passage of Fig. [3j3. 
As a result, noise in the period of oscillation of r] = 10% 
turns P(I) into a monotonously increasing function (see 



Fig. [3]3). This also explains why a small noise in the 
oscillation period drastically extends the characteristic 
time before breakdown of wave propagation in Fig. [2p- 
E. 



DISCUSSION 

The inflammatory response is initiated by the NF-^B 
regulatory network. The dynamics of NF-^B response is 
very peculiar when measured on single cell level. While 
the oscillatory dynamics of NF-^B in single cells have 
been the focus of extensive research [4j [13] , little atten- 
tion has been given to the the role of the lower secondary 
amplitude and the noise in oscillations. 

Rather than seeking the single-cell explanations, we 
have taken the perspective of spatially distributed cells, 
which collectively propagate cytokine signals through the 
tissue, thereby recruiting white blood cells from the blood 
stream to the site of infection. This approach led us to 
several interesting findings: a) The system is more sensi- 
tive to noise in refractory period than in oscillation pe- 
riod, b) Moderate noise in secondary oscillations stabi- 
lizes concentric wave patterns in the presence of noise in 
the refractory period. The effect is stronger with higher 
secondary amplitude. However, the secondary amplitude 
should never be so high that a group of cells in their sec- 
ondary oscillation can re-excite a cell to its first peak, as 
this potentially would create numerous sources of prop- 
agating waves in the system. This may help to explain 
why the experimentally observed secondary amplitude is 
indeed much smaller than the initial peak, ~ 20%. This 
prediction can be tested experimentally. The propaga- 
tion of waves could be measured in a 2D culture of mam- 
malian cells with a fluorescent reporter fused to NF-kB 
[4]. In this case nuclear NF-kB will serve as a proxy for 
cytokine production. The secondary amplitude is con- 
trolled by A20, such that higher expression of A20 re- 
sults in lower amplitude. There already exist cell lines 
where the amounts of A20 can be tuned externally [23] . 
Combining the two: fluorescent microscopy of 2D culture 
of cells with externally tunable A20, one can test our 
predictions that higher amplitude would stabilize wave 
propagation. 

c) Increasing noise in refractory period increases the 
chance of self- sustained excitations, thus decoupling the 
excitations from the original source. This situation is in 
effect very similar to chronic inflammation - inflamma- 
tion which persists and re-occurs even after the source of 
damage has been removed [24]. The system can show a 
very rich, non-monotonous behavior in how it turns in- 
sensitive with increasing refractory noise. Remarkably, 
noise in secondary oscillations postpones the onset of in- 
sensitivity, thus rendering system more robust to low 
noise in refractory period. Our findings related to the 
effect of noise in the refractory period are limited to the 
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inflammatory response alone, but are general for other 
excitable media models. 

Overall, our investigation shows that both of the ex- 
perimentally observed features: the lower amplitude of 
secondary oscillations and the noise in the period of os- 
cillations, are crucial for stabilizing wave patterns as well 
as maintaining system sensitivity. 
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